Factorization of Numbers with the temporal Talbot effect: Optical implementation by 

a sequence of shaped ultrashort pulses 
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We report on the successful operation of an analogue computer designed to factor numbers. Our 
device relies solely on the interference of classical light and brings together the field of ultrashort 
laser pulses with number theory. Indeed, the frequency component of the electric field corresponding 
to a sequence of appropriately shaped femtosecond pulses is determined by a Gauss sum which allows 
us to find the factors of a number. 
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In 1836 Henry Fox Talbot used "a magnifying glass 
of considerable power" [l[ to investigate the interference 
pattern of light emerging from a diffraction grating pro- 
duced by Joseph Fraunhofer. Talbot noticed "a curious 
effect" : The interference patterns in planes parallel to the 
grating repeated themselves periodically as the distance 
between the plane and the grating increased. Almost 
fifty years later this self-imagingeffect was rediscovered 
and explained by Lord Rayleigh 2] . Today the Talbot ef- 
fect Q manifests itself not only in electromagnetic waves 
but also in matter waves with applications ranging 
from the observation of interferences in Ceo molecules [5[ 
to lithography (6j. In the present paper we report on a 
modern day variant of the Talbot effect using appropri- 
ately shaped femtosecond laser pulses and use it to factor 
numbers. 

Factorizing numbers is an important problem in net- 
work as well as security systems 3] • Many attempts have 
been made to use quantum systems to dramatically in- 
crease the efficiency. However, even today the challenge 
remains and only small numbers Q have been factorized 
using a quantum algorithm. 

At the same time, other physics-based methods for fac- 
torizing numbers have been proposed One of them 
relies on the properties of the truncated Gauss sum 
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consisting of M + 1 terms and N is the number to be 
factored. The argument I scans through all integers be- 
tween 2 and VA for possible factors. When I is not a 
factor, the quadratic phases oscillate rapidly with m and 
the sum takes on small values. When I is a factor, then 
all the phases are multiple of 27r and the sum is equal to 
unity. 

The proposed implementations of are based on 

multipath interferences Q . Each path produces one term 
in the Gauss sum. The difficulty is to find a system which 
is experimentally accessible and in which the required 
phase in Eq. (|TJ) is obtained by a simple variation of a 
physical parameter. 



So far this strict condition has not yet been fulfilled. 
Nevertheless, several experiments in which each phase 
of the Gauss sum is separately computed have recently 
succeeded to demonstrate the ability of Gauss sums to 
factorize numbers with physical systems. In two exper- 
iments based on NMR techniques [H| the nuclear 
spins are driven by a series of radio-frequency pulses. In 
a more recent experiment, cold atoms are excited by a 
sequence of Raman 7r-pulses [l2[ • 

In the present paper we introduce an all optical ap- 
proach towards factoring numbers relying on modern 
pulse shaping technology. Indeed, the generation of ar- 
bitrarily shaped optical waveforms [HI is of great inter- 
est in a number of fields ranging from coherent control 
[l4j to information processing [2, [TH, For exam- 

ple, pulse shapers have led to an elegant implementation 
of the Grover search algorithm using Rydberg atoms as 
quantum re gist ers fl5| . Moreover, optical realizations of 
the Grover [l8l | or the Bernstein- Vazirani [l|| [2(| algo- 
rithms have been used. Our work extends this line of 
research to factoring numbers using the Talbot effect. 

Three elements determine the Talbot effect: (i) a grat- 
ing which is periodic in space and creates a periodic spa- 
tial field distribution, (ii) interference of the waves emerg- 
ing from each slit of the grating and (iii) the paraxial 
approximation of classical optics which leads to the ac- 
cumulation of quadratic phases in the time evolution of 
these waves. As a consequence, the intensity distribution 
of light on a screen is determined by a Gauss sum. 

The present implementation follows exactly this recipe 
except that it takes place in the time rather than the 
space domain. For this temporal Talbot effect [2f|, we 
consider an electric field 
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consisting of a sequence of short pulses approximated by 
delta functions. The pulses of carrier frequency lol and 
phases 9 m appear at times T m and the weight factors w m 
guarantee that the energy of the pulse remains finite. 
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of this pulse sequence defined by the Fourier transform 
follows from the interference of the Fourier components 
of the individual pulses 
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with Alu = lu — lul. 

By imprinting appropriate phases on the pulses with 
pulse shapers we can obtain the quadratic phases char- 
acteristic of the Talbot effect. For example, the special 
choice 9 m = —2"nm 2 N/l with times r m = toT and the 
weight function 
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yields for Alu — 0, that is lu — lul, the Gauss sum Eq. (JlJ. 
This Gauss sum is thus directly calculated by multipath 
optical interferences between the optical pulses. 

The laser system is a conventional Ti: Sapphire laser 
delivering pulses of tl = 30 fs at 805 nm with 80 MHz 
repetition rate. The laser pulses are shaped with a pro- 
grammable 640 pixels phase and amplitu de p ulse-shaper 
offering a shaping window of T w ~ 30 ps [22j . 

In order to generate at once the shaped pulse sequence 
required by Eq. the complex spectral mask 
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is applied with the pulse shaper to modify the Fourier 
Transform limited input laser pulse: E out (lu) = 
Hg (lu) Ei n (lu). Each term of the sum in Eq. ^ is there- 
fore produced by one ultrashort pulse delayed by r m and 
with an extra phase shift 9 m . Here we choose T — 200 fs 
in order to produce a sequence of well separated pulses. 

The interference produced by the pulse sequence is sim- 
ply analyzed with a high resolution spectrometer. We 
measure the spectral intensity at the central wavelength 
Al = 2n c/lul and thus retrieve the Gauss sum for each /. 
The experiment is performed for I ranging between 2 and 
y/~N in order to discriminate factors from non- factors. 

Figure [T] shows two typical spectra obtained for N = 
105 and I — 3 or 9 with M + 1 = 4 pulses and after 
averaging over ca 10 s. The complex structure reflects 
the multipulse interferences and underlines the requested 
high resolution (about 0.06 nm). The normalization is 
obtained from full and minimal shaper transmission. 

Several numbers have been factorized with this 
method. In Figure [2] we display the results of our op- 
tical implementation of the factorization scheme based 
on Gauss sum for N = 105 = 3*5*7 obtained with a 
four pulse sequence (a), and for N — 15251 = 101 * 151 
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FIG. 1: Spectra of a M + 1 = 4 pulses sequence (N = 105). 
I = 3 (a) and I — 9 (b). The vertical line represents lul- 



with a nine pulse sequence (b). The first example con- 
sists of the product of twin primes whereas the second 
consist of quite far primes allowing to test the validity of 
the method. 
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FIG. 2: Experimental realization of factoring using a sequence 
of shaped ultrashort pulses: (a) N = 105 = 3*5*7 with 4 
pulses, (b) N = 15251 = 101 * 151 with M + 1 = 9 pulses. 
Experiment (dots); Theory (crosses). 



The experimental data indicated by black dots are 

2 

compared with the expected values (I) depicted 

by crosses and the agreement is very good, particularly 
for the factors whose Gauss sum comes out very clearly. 
The experimental contrast is in general smaller than 
expected. This reduction could be due to several ex- 
perimental limitations : (i) Our shaper is pixellated in 
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the spectral domain and therefore introduces temporal 
replica. These replica are separated by 35 ps and are 
particularly broad and weak due to the nonlinear dis- 
persion in the mask plane [23l . l2~ij which was carefully 
calibrated. This time window of 35 ps is restricted down 
to 28 ps by the effect of the gaussian envelope due to 
the spatial beam profile in the mask plane [12, HH ■ Its 
consequences are limited here by working on only a frac- 
tion (3 to 7 ps) of the shaping window, (ii) Another 
consequence of pixellation is the hole in amplitude as- 
sociated to large phase steps between consecutive pixels 
[25| which may induce small distortions as compared to 
the ideal transmission Hg(u>). (iii) The main limitation 
to the extinction ratio (currently of 20 dB) is due to the 
gaps between pixels in the LCD (3% of the pixel width) 
adding a non programmable pulse at t = 0, which partic- 
ipates also to this loss of contrast. This contribution is 
difficult to compensate and produces undesired interfer- 
ences with the pulse train[22]. (iv) Finally the resolutions 
of both pulse shaper and spectrometer limit the ultimate 
contrast which can be achieved. Both are carefully cali- 
brated following the procedure described in [23] ■ 

A key issue in the efficiency and reliability of this 
scheme is the choice of the truncation parameter M of 
the Gauss sum. This question is closely related to the 
phenomenon of ghosts factors [2^]. Indeed, for certain 
integer arguments I, the Gauss sum can take values close 
to unity even when I is not a factor of N. Ghosts can be 
suppressed [26| below the threshold of 1/ y/2 by choosing 
M ~ 0.7-^jV. 

The example N = 19043 = p(p + 2) with p = 137 
is perfectly suited to test the predictions of Ref. [261 ] 
concerning ghost factors. In this case, N consists of the 
product of twin primes which are approximately equal 
and of the order of y/~N = 137.996. In this way we can 
test our method at the upper boundary y/N of our set 
of trial factors. For this purpose we first note that for 
any integer number N consisting of the product of twin 
primes the elementary relation N = p(p+2) = (p+1) 2 — 1 
yields the approximation p + 1 = \/N together with the 
decomposition N/(p+ X) = (p + 1) — As a 

result the truncated Gauss sum reduces to 
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Since 1 <C M and 1 <C (p + 1) we can a ppr oximate this 

sum by a Fresnel integral which yields [2 6J the scaling 

M CX y/pTT = \VN. 

In Fig. [3] we display by crosses the exact sum 
2 
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as a function of M, We note the slow de- 
cay and the oscillations due to the Fresnel integral. Solid 
dots representing our measurements follow this behavior. 
The general trend is well reproduced. However the exper- 
imental uncertainties do not allow to reproduce fully the 
expected oscillations. Moreover, we find the predicted 
threshold M ~ 0.7</l9043 ~ 8. In the insert, an exper- 
imental realization of factoring N — 19043 = 137 * 139 



with a 9 pulses sequence is shown as an example. Theory 
and experiments are also in excellent agreement. 
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FIG. 3: Supression of the ghost factor I = 138 of N = 19043 = 
137*139 for increasing number of pulses. Insert : Experimen- 
tal realization of factoring iV = 19043 with M + 1 = 9 pulses. 
Experiment (dots); Theory (crosses). 



Our work clearly demonstrates that we can use shaped 
femtosecond pulses to implement Gauss sums and fac- 
tor numbers. However, many generalizations offer them- 
selves: (i) So far we have only made use of the phases 
9 m in the frequency representation Eq. (01 of the electric 
field. The second contribution to the phase, that is the 
product r m Au> did not enter since we set Auj — 0. (ii) 
Since we have only the single parameter 9 m at our dis- 
posal, the number N to be factored and the trial factor 
I cannot be varied independently, (iii) Finally we have 
pursued a sequential rather than a parallel approach. In- 
deed, we have only used a single spectral component. 

The activation of the so far unused phase r m Acu solves 
all three problems. Since now we have two parameters 
we can encode N in r m and I in Alo. By recording the 
complete spectrum we achieve a massive parallelism. 
The choice of 8 m = and r m = 2-Km 2 Na with the nu- 
(7) merical constant a also yields the Gauss sum -4^^ and 
illustrates this new approach. Here the spacing between 
pulses increases quadratically and I is inversely propor- 
tional to Au such that a single spectrum directly contains 
all the information. 

However, some remaining difficulties need to be over- 
come: (i) The Gaussian shape of the spectral profile leads 
to ponderations in the Gauss terms which have to be 
taken into account, (ii) The variation of I between 2 and 
yJ~N i. e. on several orders of magnitudes puts severe 
constraints on the spectral resolution necessary to carry 
the experiment, (iii) Finally, the number of pulses is lim- 
ited by yfTuj/'STL ~ 10 with our present set-up. 

The quadratic spacing of the pulses required in the 
above approach might represent a severe problem. The 
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choice 9 m = and r„ 
the Gauss sum [9| 



2w(m + m 2 /N) which leads to 



W m exp 
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might be an interesting way around, since it also allows 
us to factor numbers. In contrast to the truncated Gauss 
sum which only needs to be recorded at integer ar- 

guments, the sum Sn relies on a continuous argument. 
Here we need the complete spectrum. In return Sn dis- 
plays interesting scaling properties which enable us to 
use the interference pattern for N to factor the number 
N' by rescaling the frequency axis [27l |. 

We conclude by noting that our implementation is 
closely connected to the Talbot effect in a harmonic os- 
cillator [28j . Indeed, an initial wave function consisting 
of an array of sharp maxima located at integer multi- 
ples m of a period d accumulates quadratic phases in its 
time evolution. On first sight this behavior is surprising 
since the energy spectrum of the harmonic oscillator is 
linear. However, due to the quadratic dependence of the 
energy on the position, the maxima at md translate into 
quadratic phases (md) 2 . A similar behavior was noted 
[29( in the quantum carpets woven by a wave packet mov- 
ing in a harmonic oscillator and consisting only on energy 
eigenstates with quadratic quantum numbers. 

Although the electromagnetic field represents a har- 
monic oscillator it might be easier to realize a factor- 
ization scheme based on this effect using a mechanical 



oscillator. A laser cooled atom into an optical lattice 
and prepared in its motional ground state offers a pos- 
sible realization. An absorption grating [30] produced 
by a standing light wave can prepare the periodic array 
of narrow wave packets. Moreover, the coupling of the 
center-of-mass motion to a quantized standing light field 
[3l| can introduce entanglement into the Talbot effect 
making factorization with entangled Gauss sums viable. 

In summary, we have factorized numbers through the 
implementation of a Gauss sum with optical interferences 
produced by a sequence of shaped short laser pulses. This 
work opens the route to further promising developments 
based on the wide flexibility offered by optical interfer- 
ences. 
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